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We show that anomalous diffusion arises in two different models for the motion of randomly 
forced and weakly damped particles: one is a generalisation of the Ornstein-Uhlenbeck process 
with a random force which depends on position as well as time, the other is a generalisation of 
the Chandrasekhar-Rosenbluth model of stellar dynamics, encompassing non-Coulombic potentials. 
We show that both models exhibit anomalous diffusion of position x and momentum p with the 
same exponents: {x^) ~ Cxt^ and (p^) ~ Cpt'^'"'. We are able to determine the prefactors Cx, Cp 
analytically. 
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INTRODUCTION 

In many systems the growth of a dynamical variable 
X with time t satisfies {X"^) ~ Ct" where the angu- 
lar brackets denote averaging. The process is said to be 
anomalous diffusion if a ^ 1. Anomalous diffusion may 
be a consequence of a power-law built into the dynami- 
cal process, such as in Levy flight models [1| , or it may 
be an 'emergent' property, where the anomalous expo- 
nent a is not a direct consequence of power laws which 
are built into the model. The latter case is more inter- 
esting, because non-integer exponents do not feature in 
the fundamental laws of physics, but there are relatively 
few models where 'emergent' anomalous diffusion can be 
analysed exactly. In this paper we describe two physi- 
cally natural models for the diffusion of a particle which 
is accelerated by random forces. If the damping is suffi- 
ciently weak the particle can exhibit anomalous diffusion, 
having universal exponents, the same for each model. 

Our two models are generalisations of two classic mod- 
els for diffusion processes. The first is an extension of 
the Ornstein-Uhlenbeck process [2], in which a particle 
is subjected to a rapidly fiuctuating random force, and 
is damped by viscous drag. In the generalised Ornstein- 
Uhlenbeck process the random force depends upon the 
position of the particle and time and is derived from a 
potential $(a;,i). Earlier works analysed this model in 
detail for one spatial dimension [3|, |j| . This model ex- 
hibits anomalous diffusion in one dimension. Here we 
discuss higher spatial dimensions, where the mechanism 
for anomalous diffusion is significantly different, as was 
suggested (for a closely related model) in j5|, Q. Here 
we obtain exact formulae for the momentum distribution 
of the generalised Ornstein-Uhlenbeck process in two and 
three dimensions (for a particle which is initially at rest). 
We use these to obtain precise asymptotic formulae for 
the growth of the second moment of the coordinate: the 



second moments scale as (p^) ^ t^/^ and {x^) ^ t^ in the 
anomalous diffusion regime (the exponents are the same 
as those obtained in [6| ; we obtain the prefactor exactly) . 
We also discuss an extension of the Chandrasekhar- 
Rosenbluth model for diffusion [7|, |8|, in which a test 
particle interacts with a gas of point masses via a pair 
potential. The interaction should cause small changes of 
momentum, which can be modelled as a diffusion process. 
Usually the interaction is gravitational, and the applica- 
tion is to the motion of stars in galaxies, but here we 
simplify the problem by considering a non-singular weak 
interaction potential. We show that, surprisingly, the 
diffusion tensor has the same form as for the generalised 
Ornstein-Uhlenbeck process, and that consequently there 
is anomalous diffusion with the same universal exponents. 
For this model too, we derive diffusion coefficients for this 
model precisely in terms of the microscopic parameters. 

The anomalous diffusion effect which we describe is 
analysed by introducing a diffusion process describing the 
fiuctuations of the momentum p of a particle in response 
to a spatially and temporally fluctuating random poten- 
tial. The diffusion coefficient of this process, D{p), is a 
function of the momentum of the particle. We remark 
that this approach to formulating the equation of motion 
was first introduced by Sturrock ^] , and that similar de- 
velopments appeared later in the mathematical literature 
(see [1Q!| and references cited therein). Our paper is the 
first work to give a solution to the generalised Ornstein- 
Uhlenbeck process in two or three dimensions. 



GENERALISED ORNSTEIN-UHLENBECK 
MODEL 



We consider a particle of mass m with momentum p 
subjected to the generalised Ornstein-Uhlenbeck process 



I . This is described by equations of motion 
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P= -1P + f{x,t) 



(1) 



where x is the particle's position. The particle experi- 
ences two types of forces: a drag force — 7P with 7 being 
a damping rate and a random force f{x,t). Unlike the 
classic Ornstein-Uhlenbeck process, where the force de- 
pends only upon time, in our generalised model the ran- 
dom force depends upon the position as well. We assume 
that f{x,t) is a force derived from a random potential 
varying in time and space, i.e. f{x,t) = — V$(a;,i), 
where ^{x,t) has statistics 

($(x,t)) = , mx,m{x',t')) = C{\x ~ x'\, \t - t'\) . 

(2) 
The correlation function C{x, t) has temporal and spatial 
scales, r and ^, respectively. We consider the case where 
7T ^ 1, implying that the momentum satisfies a diffu- 
sion equation. We define a momentum scale pq — m^/r, 
such that if |p| 3> po the force experienced by the parti- 
cle decorrelates much more rapidly than the force expe- 
rienced by a stationary particle. We consider the limit 
where \p\ 3> Po, which is realised for weak damping. The 
dynamics of ^ can be described by a diffusion equation 
for the probability density of the momentum, P{p, t): we 
now consider how to derive this diffusion equation. 

The dynamics of the momentum can be approximated 
by a Langevin process. Small increments of components 
Pi of the momentum vector p may be written as 



Spi = -IPiSt + Swi 



(3) 



where Swi is the impulse exerted by the i-th component 
of the force / on the particle in time St, 

j-to+St 

Swiito) = \ dii fi{x{ti),ti) 

rto+st 

dii Mpti/m,h) + 0{St^) . (4) 

The Langevin process ([3]) is equivalent to the Fokker- 
Planck equation describing time evolution of the proba- 
bility density of momentum P{p, t). In order to construct 
the equation we need to know drift and diffusion coeffi- 
cients, Vi — {6pi)/6t and Dij — {Spidpj) /2St respectively. 
Using the definition of the increment dwi we find 

rSt pSt 

{5wi6wj) ^ j dti dt2{fi{pti/m,t)fj{pt2/m,t2)) 
Jo Jo 
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~ 6t / dh (/,(O,0)/j(pti/m,ii)) 
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(5) 



(the second step is justified when 5t is large compared to 
T but small compared to 7""'^). 



The components Dij of the momentum diffusion tensor 
D depend upon the direction of the momentum. For the 
case where the momentum is aligned with the x-axis (that 
is, where p = pei, the coefficients of the diffusion matrix 
are expressed in terms of the correlation function of the 
potential as follows: 
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(In three dimensions D^z — Dyy and D^z — Dyz = 0.) 
For other directions of the momentum, the elements of 
the momentum diffusion tensor can be obtained by ap- 
plying a rotation matrix. If O is a rotation matrix which 
rotates the momentum vector p into p' = Op = pei 
which is aligned with the x-axis, then the elements D[, 
of the diffusion matrix in the transformed coordinate sys- 
tem are given by ^ . The diffusion matrix in the original 
coordinate system is D = O D' O""". 

Having considered the diffusive fluctuations of the mo- 
mentum, we now consider its drift, (Spi) — ~^piSt + 
(dwi). Expanding fi(x,t) we obtain (in d dimensions) 



Mx,t)^Mo,t) + 



^ dx, " 



,(i) (7) 



where Xj{t) can be written as a solution of ([T]), 



1 



it) = - dti / dt2 exp[-7(ii - ts)] 



X fj{pt2/m,t2) . (8) 

Combining together Q, ^ and ^ we obtain 
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dh / dt2 / dig exp[-7(i2 - h)] 
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Ij{ph/m,h) 
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We approximate exp[— 7(^2 — ts)] by unity for 7T <C 1, 
and use the assumption that (5i 3> t to obtain 
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Taking account of the fact that the drift coefficients con- 
tain derivatives of the diffusion coefficients, we obtain the 
Fokker-Planck equation, 

dP d f d\ 

_._. -,p + D(p,- P. (11) 



We are primarily interested in the case of weak damp- 
ing, where the typical momentum satisfies \p\ ^ Po- In 
this limit, the diffusion coefficients have an algebraic de- 
pendence upon p = \p\. For diffusion of the momentum 
vector parallel to its direction, when p ^ po the leading 
term is of order p^^ and the result is given by 
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For diffusion of p in a direction perpendicular to its di- 
rection, we have 
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2po 
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(13) 

Note that when p ^ po, the diffusion of the direction 
of the momentum vector is much more rapid than diffu- 
sion of its magnitude. This makes the behaviour of the 
generalised Ornstein-Uhlenbeck process in two or more 
dimensions very different from its behaviour in one di- 
mension. 

Results equivalent to (fT2|) and (TTSl) were obtained in 
[a, |6| in an analysis of a closely related model. 



GENERALISED 
CHANDRASEKHAR-ROSENBLUTH MODEL 

In this section we consider motion of a particle travel- 
ling through an infinite homogeneous population of back- 
ground particles. We assume that the test particle inter- 
acts with the background particles, and the background 
particles do not interact with each other. As the test par- 
ticle moves, its interaction with each of the background 
particles causes small changes of its velocity. When the 
number of the background particles is very large, the ve- 
locity of the test particle changes rapidly and in an un- 
predictable way, so that its motion can be described by a 
diffusion process. The test particle could be a star mov- 
ing in a galaxy interacting with the background stars (the 
interaction between the background stars is not consid- 
ered) , so that the force of interaction is gravitational and 
thus proportional to the inverse square of the distance 
r between stars. This problem was originally studied by 
Chandrasekhar [7| , who found that the test particle expe- 
riences a gradual decrease of the velocity in the direction 
of motion. This phenomenon is called 'dynamical fric- 
tion'. For 'slow' particles (with velocities much smaller 
than some representative velocity scale) this deceleration 
is proportional to the velocity of the particle v (analo- 
gous to the Stokes's law for a drag force for a particle 
is a viscous medium). For sufficiently 'fast' particles the 
deceleration is proportional to v~'^. 

Subsequently, Rosenbluth et al. [8| studied the dif- 
fusion of momentum of the test particle in this model 
in greater depth. They found that in the 'fast' regime 



the diffusion coefficient of the momentum in the direc- 
tion parallel to the direction of motion is proportional 
to v^^ , while the diffusion coefficients in the plane per- 
pendicular to the direction of motion is proportional to 
v~^. These dependences are equivalent to the momen- 
tum dependences of the radial and transverse fluctuations 
of the momentum in the generalised Ornstein-Uhlenbeck 
model, obtained in equations (|12|) and (fTB]). The ex- 
pressions for these diffusion coefficients obtained in [8| 
contain logarithmic terms due to the long-ranged na- 
ture of the Coulombic potential, which makes it difficult 
to write down precise formulae. In order to illuminate 
the relation between the generalised Ornstein-Uhlenbeck 
model and the Chandrasekhar-Rosenbluth model in the 
simplest context, in this section we consider the latter 
model for the case when the interaction is described by 
a short-range potential U{r) of some rather general radi- 
ally symmetric form. We obtain precise expressions for 
these diffusion coefficients and show that the results for 
the scaling of the diffusion coefficients are the same as in 
the generalised Ornstein-Uhlenbeck process. This leads 
to the same anomalous diffusion behaviour in both mod- 
els. 

We shall only discuss the two-dimensional case for sim- 
plicity and proceed as follows (some of the presentation 
adapts the discussion of the Chandrasekhar model in 
HJ). We first calculate the change of the velocity of a 
test particle due to the encounter with a stationary back- 
ground particle. We denote components of this change 
by Avn and Av'j_, in the directions parallel and perpen- 
dicular to the initial velocity of the test particle. Next, 
we consider the change of the velocity of the test parti- 
cle for the case when the background particle propagates 
with velocity vi,. We denote components of this change 
by Awii and Avj_. Using geometrical arguments, we then 
express Av± and Af || in terms of Aii^ and Awf, . 

Let us consider a test particle of mass m moving in the 
horizontal direction with the initial velocity Vq = (Vb, 0). 
The test particle interacts with a background particle of 
mass Af , which is initially at rest. After an encounter 
the test particle propagates with velocity Vi described 
by its magnitude Vi and the polar angle ^i , and the back- 
ground particle moves with velocity V2 described by its 
magnitude V2 and the polar angle ^2 (see Fig. [1]). 

The changes of the velocity of the test particle in the 
direction parallel and perpendicular to Vq are 



AW|j = ViCOS^i -Vq, 

Aw^ = Vi sin ^1 . 



(14) 



The conservation of momentum before and after the en- 
counter yields 



mVo = mVi cos ^1 + MV2 cos ^2 
= mVi sin ^1 -t- MV2 sin ^2 



(15) 
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FIG. 1: The encounter between a test particle (black) and 
a background particle (shaded, red online). The interaction 
causes a small change of the velocity of both particles. 



and from the conservation of energy we obtain 



(16) 



This enables us to write an equation for Vi: 

V^(m + M) ~ 2mVoVi cos^i + V^{m - M) = . (17) 

We assume that the encounter induces only a small 
change of the direction of motion of the test particle, 
so that ^1 can be taken as being small. In this approxi- 
mation the solution of equation (|17p is 



Vi ^ Vo{l - aCl) + 0{Ct) 



(18) 



where a = m/{2M). We substitute Vi from (|18p into 
equation (fTl]) and obtain 



Awf, 



Ai 






(19) 



where [3 = a+\/2 = {m+M)/2M. In the small-angle ap- 
proximation ^1 is defined by the change of the momentum 
in the perpendicular direction, so that ^i ^ Ap±/{mVo)- 
The change of the momentum Ap± is determined by the 
force of interaction between particles separated by dis- 
tance r with magnitude /(r) = —dU{r)/dr, so that Ap± 
can be written as 



ApA 



dt^[r{t)]smx{t) (20) 

dr 



where r{t) is a distance between the particles and x(i) 
is an angle between Vq aiid the vector connecting the 
particles. We define an impact parameter b to be the 
initial distance between the test and background particles 
along the axis perpendicular to Vo- From figure[l]we find 
that r(i) = Vx{t)^ + 6^ and sinx(t) = b/^/x{tf + P, 
where x(t) is a coordinate of the test particle along the 
direction parallel to Vq. Changing the variable x{t) — 
Vot, in the weak-scattering limit where the defiection is 
small we obtain 



1 f°° , dU, r 

Ap± sa -— / dx -5- V2 
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(21) 



If we denote an integral 

we obtain 

Using this relation we find 



mV^ 



(22) 



(23) 



Av', 



Av\ = 



mVo 



(24) 



Thus, the contribution to the change of the velocity of the 
test particle due to a single encounter is proportional to 
Vq~ and Vq~ in the directions parallel and perpendicular 
to Vq , respectively. Averaging over collisions with many 
particles is expected to add another factor of Vq, as the 
particle propagates with this velocity (see also the deriva- 
tion below). This suggests that the second moments of 
the change of the velocity scale as Vq~^ and Vq"^ in the 
directions parallel and perpendicular to Vq, respectively. 
While the latter result is consistent with the behaviour 
of the diffusion coefficient in the generalised Ornstein- 
Uhlenbeck process, the former result is different, as in 
the previous model the result is (Av?,) ~ v~^. However, 
if the background particles are not stationary, these esti- 
mates must be corrected, as shown below. 

We assume that the background particle moves with 
velocity Vh, in which case the discussion above is valid if 
Vo is a relative velocity of the test particle in the frame 
of reference moving with the background particle. The 
velocity of the test particle in a fixed frame of reference is 
therefore Vq = Vq + vi,- We are interested in the changes 
of the velocity of the test particle in the directions parallel 
and perpendicular to Vq. We denote these Awy and Av± 
and deduce from figure [2] 



Aw|| = Aw|| cosfi -)- Aw^ sinJl, 
Av± = Av\ cos il — Au[| sin fl 



(25) 



where 17 is an angle between Vq and Vq. These relations 
are equivalent to the rotation of the coordinate system 
by angle fl. Assuming that vq^ Vb, where Vb = \vb\ and 
''^o = \vo\, we obtain 



r^- 



t'fc A Vq ^ Vb±_ 
V^ ~ Vq 



(26) 



and Vq ^ Vq. In this approximation, to the leading order 
in Vq, we have 



Awii 



Av±_ 



I{b)vbA 



m 

mvQ 



(27) 




FIG. 2: Geometrical construction illustrating the changes of 
the velocity of the test particle Av parallel and perpendicular 
to the relative velocity Vo {Avi and Av'^, solid lines, blue 
online) and the velocity of the test particle vq (Aii|| and Av±, 
solid lines, red online). 



We now imagine that the test particle is travelling 
through an infinite homogenous population of the back- 
ground particles with the spatial density number n (mea- 
suring a number of particles per unit area) and the prob- 
ability density of the velocity is f{vi,). Let diV be the 
number of background particles it encounters in time At 
with velocity Vf, in a volume element of velocity space dvj, 
and impact parameter between b and b + db. This is the 
number of particles in two thin stripes, each of width db 
and length equal to the distance travelled by the particle 
in At, multiply by the probability /(f b) dv/,. We have 



dN - 2nVaAtdb x f{vb)dvb 



(28) 



In order to obtain the total contribution of many back- 
ground particles with different impact parameters and 
velocities, we integrate over b and v^, 
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2Af 
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m'^v^ Jo 



dbl\b) dVb f{Vi,)\Vb^\' 
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(29) 



Here we used the assumption that U{r) is a short- 
ranged, allowing us to let the upper limit of the inte- 
gral over b approach infinity. In the case of the origi- 
nal Chandrasekhar-Rosenbluth model, an upper limit to 
the impact parameter must be introduced because of the 
long-range interaction between theparticles. This leads 
to logarithmic correction terms (7|, |8| . 

Equations (f^. describing the velocity increments for 
the Chandrasekhar-Rosenbluth model has the same scal- 
ing (as a function of vq) as scaling of the diffusion co- 
efficients for the generalised Ornstein-Uhlenbeck model 
(as a function of p; see equations (fT2|) and (fTS])). This 
indicates that the anomalous diffusion behaviour of these 
models is equivalent. 



PROBABILITY DENSITY FUNCTION AND 

MOMENTS OF THE MOMENTUM 

Now we return to the generalised Ornstein-Uhlenbeck 
process and obtain the closed-form solution of the 
Fokker-Planck equation for a particular choice of the ini- 
tial conditions. We use this solution to obtain an exact 
expression for the growth of the moments of the momen- 
tum. 

We first consider the two-dimensional case. The prob- 
ability density for the momentum satisfies equation (jlip , 
with the diffusion coefficients given by (fT2)) and (fT3)) . We 
transform to polar coordinates and seek a probability 
density P{p,9,t), and consider the case when the par- 
ticle is initially at rest, so that the initial condition is 
P(j),9,0) — S{p). This circularly symmetric solution, 
P = p{p,t), satisfies 

dp 
dt 
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(30) 



By analogy with the solution of the one-dimensional gen- 
eralised Ornstein-Uhlenbeck model, we find the following 
normalised closed- form solution of (15(11) : 



P(P, t) 
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X exp 
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(31) 



In the long-time limit the density is non-Maxwellian 
given by 



Po(p) 



57 
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-exp 
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5D3pi 



(32) 



T{2/5){5Dspl)y^ 

Using the probability density (|3T|) we determine the l-th 
moment of p, 



{P\t)) 



dpp'+'p{p,t) 



/5i^3p3y/^r[(2 4-0/5] 
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-57ty/5(33) 



We remark that an additional factor of p in the expression 
above appears as a weight in the transformation to polar 
coordinates. 

In the three-dimensional case we find a similar solu- 
tion to equation (fTO)) in the case where the particles are 
initially stationary. We write this equation in spheri- 
cal polar coordinates, and seek a spherically symmetric 
solution, P{p, 9, (f), t) = p{p, t). The solution of the corre- 
sponding equation for p{p, t) is obtained similarly to the 
two-dimensional case and we have 

v3/5 
P{P. t) 
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r(3/5) [5D3Pg(l - e-57t)]3/5 
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5i?3P3(l - e-5Tt) 



(34) 
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FIG. 3: Shows results of the numerical simulation of equa- 
tion ll]) for the motion in the two-dimensional potential 
force field. Panel a shows stationary non-Maxwellian den- 
sity function (|32p (solid line) and data from the numeri- 
cal simulation (circles). Panel b shows the second moment 
of the momentum (I33p with I = 2 (solid line) and data 
from the numerical simulation (circles). The dashed line 
shows the slope t^'^ and dotted line indicates time 7~^ at 
which the density becomes stationary. The results are for 
the case of a Gaussian correlation function of the potential 
C{x,t) = cr^exp[-xV(277^) - tV(2T^)] with a = 15, r = 0.1 
and ri — 0.1. The other parameters were m — 1 and 7 = 0.01. 



This determines moments of the momentum 

1/5 
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(35) 



For both two- and three-dimensional cases we obtain that 
at short times the variance of the momentum grows as 



(p2(i))^i2/5 



(36) 



Thus, at short times the momentum diffuses anomalously 
with the same exponent as in the one-dimensional model 
[a, y . The results for the stationary probability density 
and diffusion of the momentum in the two-dimensional 
case were verified by a numerical simulation, documented 
in figure [31 



SPATIAL DIFFUSION 

In this section we find the mean-square value of the 
displacement of a particle which starts at the origin: 
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where 6 is an angle between p(ii) and p{t2)- We recall 
that when the force is the gradient of the potential, we 
have Dxx ^^ Dyy for p ^ po implying that the corre- 
lation of the angle vanishes much more rapidly than the 



correlation of the magnitude of the momentum. We can, 
therefore, perform the averaging in (1371) by first integrat- 
ing over the correlation function of the angular variable, 
with the momentum held fixed, and then finally perform- 
ing the averaging over fiuctuations of the momentum. 

In the two-dimensional case the probability density 
P{9, t) of 6 satisfies the diffusion equation on a circle 
with the initial condition P{9, 0) — S{d). The solution is 
Gaussian: 
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(38) 



where V = Dipo/p^{ti). Using this probability density 
we calculate the expectation value 

(cos6') =exp(-X'|i2 -til). (39) 

We thus obtain from equations (1571) and (15^1) 
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Introducing a new variable T ~ ti — t2 wc have 
2 



(i^(or) 



™^ JO 



dti / dT 



j^ dp pip,h)p^ih)cxp (-^r 



/ dii / dpp(p,ii)/(ii) 
m^Dipo Jo Jo 



When the forcing is strong we have Dipot ^ p'^ for t ^ t, 
and therefore 
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Using equation (1331) we obtain 
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In the three-dimensional case the probability density 
of cos 9 can be found by considering the diffusion equa- 
tion on a spherical surface, with polar coordinates {9, (p), 
starting from the pole, 9 — Q. The solution P{9,(j),t) of 
the diffusion equation may be expressed as a linear com- 
bination of spherical harmonics. Because the problem 




FIG. 4: Shows results for the spatial diffusion in the two- 
dimensional potential force-field. The results from the nu- 
merical simulation (circles) are compared with equation (|43p 
(solid line). Dashed lines show the slopes t^ and t and dotted 
line indicates the time 7"^^. The parameters of the simulation 
are the same as in figure |3] 



has a rotational symmetry, the solution is independent 
of the azimuthal angle 0, and it may be written as: 



P(0,0,i)=^Aiexp 



(=0 



p3 



P; (cos 61) (44) 



where Pi{z) is a Legendre polynomial of degree I. Using 
the orthogonality relations for Legendre polynomials, in 
view of the initial condition cos 6 = 1, we obtain Ai = 
{21 + l)/2. Also, the quantity that we wish to average 
is itself a spherical harmonic: cos6' = Pi (cos 6'), so that 
only the I = 1 term in (|44p contributes to the correlation 
function. Hence we obtain 



ballistically. The results are consistent with a short-time 
asymptotic behaviour of the undamped particle obtained 
in [6[ for d > 1 . The long-time behaviour is naturally dif- 
fusive, (|a;(t)p) ^ t. In Fig. |4]we show the comparison 
of the analytical and numerical results for (|a;(t)p) for 
the case of motion in the two-dimensional potential force 
field, illustrating the short-time ballistic diffusion. 



SUMMARY 

We have investigated generalizations of two classical 
models for diffusion of a particle accelerated by random 
forces. We discussed a generalization of the classical 
Ornstein-Uhlenbeck process where the force depends on 
the position of the particle as well as time. We also mod- 
ified the Chandrasekhar-Rosenbluth model by consider- 
ing motion due to a short-range interaction potential. 
Although both models are described by different micro- 
scopic equations of motion, surprisingly, they have the 
same scaling of the diffusion coefhcients, leading to the 
same short-time asymptotic dynamics. 

We solved the Fokker-Planck equation for the gener- 
alised Ornstein-Uhlenbeck process exactly in two and 
three dimensions, building upon our earlier analysis of 
the one-dimensional case in [3|, Ij] . We have shown that 
this dynamics is characterised by anomalous diffusion 
of the momentum, with the variance which scales as 
(p^) ^ t^/^. At long time, the distribution of the momen- 
tum has been found to be non-Maxwellian. The second 
moment of the displacement grows ballistically at short 
times, that is (a;^) '-^ t^, in accord with a surmise made 
by Rosenbluth for a closely related model [6[ , and at long 
time a simple diffusive behaviour of the displacement is 
recovered. 



{cos0{\t2 - ill)) = cxp(-3P|t2 - til) , (45) 

where V = poDi/p^ is the same as in the two-dimensional 
case. Using this angular correlation function we calculate 



{\x{tr 



1 



dii / dt2 



dp pip,ti)p'^{ti) 



X exp {-3V\t2 - h\) . 



(46) 



The evaluation of the integral using the probability den- 
sity (|35|) yields 

(I^WP> = 7^^(57t + c-^-*-l). (47) 

We find that in two- and three-dimensional cases 
(|a;(t)p) ^ t'^ at short times, so that the particle diffuses 
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